APS/123-QED 



Mori-Zwanzig projection formalism: from linear to nonlinear 

Jianhua Xin£| 

Department of biological sciences, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061 

(Dated: April 17, 2009) 

The Mori-Zwanzig projection formalism is widely used in studying systems with many degrees 
of freedom. We used a system-bath Hamiltonian system to show that the Mori's and Zwanzig's 
projection procedures are mutual limiting cases of each other depending on the size of the projected 
Hilbert space. We also derived the dynamic equations of collective coordinates of a Hamiltonian 
system. 
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I. INTRODUCTION 

It is common to study dynamics of a system with a 
large number of degrees of freedom in almost every sci- 
entific field. In general it is impractical, and often unnec- 
essary, to track all the dynamic information of the whole 
system. A common practice is projecting the dynam- 
ics of the whole system into that of a smaller subsystem 
through information contraction. The procedure leads 
to the celebrated Langevin and generalized Langevin dy- 
namics. The Mori-Zwanzig formalism is a formal proce- 
dure of projection, especially for Hamiltonian systems. 
The original Mori procedure projects dynamics of the 
whole system into a sub-Hilbert space. In the literature 
application of the Mori procedure generally results in a 
generalized Langevin equation(GLE) that is "inherently 
linear in the system variables" 1]. The projection proce- 
dure developed by Zwanzig works on an enlarged Hilbert 
space, on the other hand, can lead to nonlinear GLEs 
@i S 0]- Recently Lange and GrubmuUer tried to de- 
rive the dynamic equations of some collective coordinates 
with the Zwanzig projection procedure 0. Chorin and 
coworkers discussed generalizing the projection technique 
to non-Hamiltonian systems [6]. In this communication 
we will show that the Mori's and Zwanzig's procedures 
are mutually limiting case to each other. We will also 
rederive the dynamic equations of collective coordinates, 
which differ slightly from that obtained by Lange and 
GrubmuUer. 



II. THEORY 

First we will summarize the Mori procedure. We will 
follow the discussions given in [l[ with some modifica- 
tions, and focus on the Hamiltonian systems. 

Consider a system described by the Hamiltonian, 
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i/(x,p) = ^| + y(x) 

1=1 



(1) 
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where x and p are position and conjugate momentum 
vectors. We will use mass-weighted coordinates through- 
out this paper. 

The Liouville operator L is defined as. 



^ \dpi dx. 



dxi dpi 



(2) 



For an arbitrary dynamic variable A, in this case the 
projection operator is defined as, 



PA = 



(3) 



{0(x, p)} composes the basis set for the projected sub- 
space. The inner product for two arbitrary variables A 
and B is defined as, 



{A,B) = <A^B> 

J A^B cxp {-f3H) dxdp 
J exp {—(3H) dxdp 



(4) 



where f means taking transpose and complex conjugate. 
Any dynamic variable within the subspace can be ex- 
pressed as a linear combination of the basis functions. 
The projected equations of an arbitrary dynamic vari- 
able A, which is defined within the projected subspace, 
are given by, 



where 



PLA{t) 



dsK{s) ■ (t){yi{t - s), p(t - s) + F{t){b) 



F{t) = exp(i(l - P)L)(1 - P)L^ (6) 
K(t) = 

- (F(i),L(/.)-(0, (7) 

The last equation leads to the generalized fluctuation- 
dissipation relation, and we have used the anti-Hermitian 
property of the Liouville operator. 

In practice the basis sets are usually chosen as portion 
of the coordinate vector x and the corresponding conju- 
gate momentum vector p. Then the Mori projection pro- 
cedure results in a GLE that is linear to the coordinates 
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and momentum. However, in principle this restriction is 
unnecessay. One can expand the Hilbert space to include 
high order combinations of the coordinates and momen- 
tum, as shown by the example discussed below. Here we 
examine the extreme limit of including all the possible 
Hilbert functions composed by the coordinate and veloc- 
ity (or momentum) in reduced dimension. The following 
procedure is analogous to what adopted by Zwanzig Q. 
For simplicity let's focus on a one-dimensional coordi- 
nate and its conjugate momentum, while generalization 
to higher dimensions is straightforward, 

c = /(x) = /(O) + V,/(0) • X + i Vxx./(0) : XX + . .(8) 



Vx/ • i = Vx/(0) • X + Vxx/(0) : XX 



(9) 



Therefore c and c are vectors in the full Hilbert space. 
Let's consider the sub-Hilbert space, which may still have 
infinite dimension, supported by all the possible multi- 
plicative combinations of c and c , such as (?^c(? .... 
A key observation is that these basis functions, denoted 
{0(c, c)}, compose a complete basis set for the subspace 
so any arbitrary function of (c, c) , can be fully expressed 
by the basis set. That is, for an arbitrary function 
5(x,p), 



E 



; exp(— /3H)dxdp 



= / 5exp(-/3//)(5(/ - c)5(Vx/ • p - c)dxdi^lO) 

P(c, c) J 



where 



p(c, c) = y exp(-/3i7)^(/ - c)^(Vx/ ■ p - c)dxdp (11) 

The above expression may be more familiar if the Dirac 
bra and ket notations are used. Then for the projected 
equations. 



PL-c 
PL-c 



c 

-keT 



1 d 



j exp(-/3iJ) 



(12) 



p{c, c) dc 

Vx/|P<5(/ - c)(5(Vx/ • p - c)dxdp (13) 



To derive the above expression, we performed integration 
by parts, and used the relations, 

Vx<5(c - /) = V^fdf6{x - /) = Vx/9c5(x - /) 
Vx(5(c - Vx/ • p) = Vx(Vx/-p)a£j(c- Vx/-p) 
Vp(5(c - Vx/ • p) = VxM(5(c - Vx/ ■ p) (14) 

Compared to Eq. 1131 the result derived by Lange and 
GrubmuUer has an extra term ||Vx/|P in the expression 
of p(c, c) ■ The discrepancy may come from the fact 
that the projection operator defined in [5| does not rig- 
orously satisfy P^ = P. It remains to be examined on 
how this extra term may affect the dynamics. In the case 
/ is a linear combination of x, and is chosen to satisfy 
||Vx/|P = 1 , Eq. [T3l gives the familiar relation to the 
gradient of potential of mean force. 



III. EXAMPLE 



Here we consider a system-bath Hamiltonian, 



H 



2^ 



+ 9, 



l3 



dx{t) 

dt 
dp{t) 

dt 



Pit) 

-x{t) - bx{tf 



KN{t) = ^^cos{ujjt) 

FN(t) 



E sincjji 



with the fluctuation-dissipation relation, 

< FN{t)FN{t') >o - keTKNit - t'). 



dx{t) 

dt 
dpjt) 

dt 



vit) 

-cj^xit) 



< FL{t)FL{t') > = keTKLit-t') 



>15) 



Zwanzig discussed a nonlinear GLE for the system coor- 
dinates {x,p} obtained by directly solving the equations 
of motion 1, 7|] , 



dsKN{s)p{t - s) + FN{t) (16) 



The memory kernel and the random force terms are given 
by, 



(17) 



+ E^J (^jW - 5^^(0)^ cosa;,t (18) 



(19) 



The average is over an equilibrium heat bath with the 
system constrained at {a;(0),p(0)}. By projecting to the 
Hilbert space (x, v) with Mori's procedure, one can also 
obtain a linearized GLE Jj], 



dsKL{s)p{t - s) + FL{t){2Q) 



Where Wq — ksT/ < x^ >, and the the random force and 
memory kernel terms are also related by the fluctuation- 
dissipation relation 



(21) 



However in this case, the average is over the uncon- 
strained thermal equilibrium distribution. Effects of the 
nonlinear term —bx{t)^ are contained in the renormalized 
coefficent ujq, the memory kernel, and the random force 
terms. 

In the following discussions, we will generalize the 
projection procedure of Mori by choosing a basis set 
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{x, x^,v}. Functions with even powers of x makes no con- 
tribution to the projection {{Lp, x^") = 0, n = 1,2,...). 
Therefore the lowest nolinear basis function is x^. 
First, 



Lx — p, Lx^ 



Lp = —X — hx^ — 



7^ I 2^ 



- q. 



(22) 



However, one also has. 



< Lp, > 



< X 

E 



Tl+l 



> < 



r.n+3 



> 



'T'i 



= - < X 



n+l 



> -b< 



> 



(28) 



let's calculate the normalization matrix, 



A- 



< x^ > < x"^ > < xp> 

< x'^ > < x^ > < x^p > 

< xp > < x^p > < p'^ > 

<x^ > /h - <x^ > /h 
- < x^ > /h < x^ > /h 




(23) 








24) 



Where h =< >< x^ > 



< x^ The memory 

function and the random force in the equation of motion 
of x vanish, which can be seen from, 

Lx = {{Lx,x) {Lx,x^) {Lx,p))-A^^ 

= P (25) 



One has, 
(Lp,x") 



1 



/ exp{—PH)dx 



J exp{—f3H)dx 
f exp{-(3H)da 



x"—exp{-PH)dx 
exp{~l3H)dx 



dx 

nx^^^ exp{-PH)dx 



> 



{Lp, x) = —ksT = —uJq < x^ > 



(26) 
(27) 



Therefore, 



< x^ > 



< > 



1 



(c.2 - 1) < x2 > 



-Wg < X- >< X" > 

-^(^0 - 1) < 2:" > 



Then, 



PLp{t) = -x- bx^ 



(29) 



(30) 



(31) 



One can easily show that the random force (through 
dF/dt = (1 — P)LF) and memory kernel (through Eq. 
[7]) terms are the same as those given in Eqs. [171 and [TSl . 
although in general here the average perform in Eq. [7] is 
over the unconstrained thermal equilibrium distribution. 
Therefore with the Mori projection procedure we recover 
Eqs. [TBI [T71 [T51 obtained by exact integration. Following 
similar procedure, one can show that further expanding 
the basis functions to include higher orders of x" does not 
change the projected equation form. The above results 
can also be obtained by applying Eqs. [5l I12I13I directly. 
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